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Abstract 



The long-time evolution of a system in interaction with an external 
environment is usually described by a family of linear maps 7t , generated 
by master equations of Block-Redfield type. These maps are in general 
non-positive; a widely adopted cure for this physical inconsistency is to 
restrict the domain of definition of the dynamical maps to those states 
for which 7 t remains positive. We show that this prescription has to be 
modified when two systems are immersed in the same environment and 
evolve with the factorized dynamics 7t (g) 7* starting from an entangled 
initial state. 
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1. INTRODUCTION 



The dynamics of systems immersed in large, external environments can be described 
in terms of master equations: they generate the finite time-evolution for the reduced 
density matrix, obtained by tracing over the environment degrees of freedom. [1-9] Their 
explicit form is in general rather complex, involving non-linearities and memory effects. 
Nevertheless, when the coupling between subsystems and environment is sufficiently weak 
and for times much longer than the characteristic correlation-time in the environment, 
suitable limiting master equations in Markovian form can be derived. 

These derivations are often based on ad hoc approximations, lacking mathematical 
rigor, while the final result is justified on the basis of physical considerations. Despite 
these heuristic treatments, Markovian master equations have been applied to model various 
effects in open system dynamics, ranging from quantum optics to quantum chemistry. 

It has been pointed out long ago that the heuristic derivations of the Markovian limit 
of master equations could lead in general to physical inconsistencies. [10] In particular, the 
resulting finite time-evolution described by such equations would not in general preserve 
the positivity of the reduced density matrix, with some remarkable exceptions, based on 
rigorous mathematical treatments. [1-4, 11-19] 

Although acknowledged in most subsequent literature on the subject, these incon- 
sistencies were either dismissed as irrelevant for all practical purposes [20] or cured by 
adopting further ad hoc prescriptions. [21-23] In the latter case, the general attitude is to 
restrict the action of the dynamical maps generated by the non-positive master equations to 
a subset of all possible initial reduced density matrices, those for which the time-evolution 
remains positive. This is equivalent to a suitable selection of the initial conditions for 
the starting state of the subsystem, a procedure sometimes referred to as "slippage of the 
initial conditions". On physical grounds, this effect is viewed as the consequence of the 
short-time correlations in the environment, that have not been properly taken into account 
in the derivation of the Markovian limit of the original master equation. ^ 

In the following, we shall reexamine this widely used prescription to cure possible 
inconsistencies produced by non-positive, Markovian master equations, and point out fur- 
ther potential problems of this approach. We shall deal with two identical, non-interacting 
subsystems immersed in a same environment, both evolving, in the Markovian limit, with 
the same non-positive master equation. We shall explicitly show that redefining the ini- 
tial conditions to make positive the single system time-evolution is not enough to cure all 
possible inconsistencies of the two-system dynamics. These show up when the two-system 
state that emerges after the transient due to the short-time correlations in the environment 
is entangled; therefore, in order to have a physically acceptable time evolution for the two 
subsystems when entanglement is the most likely consequence of the initial transient phase, 
the above mentioned procedure of restricting initial conditions should take into account 
also correlated states. 



' Let us point out that, instead of restricting the possible initial states, one can alter- 
natively "smooth" the initial conditions on which the non-positive dynamical map acts 
[24] ; the resulting effective map turns out to be positive. Unless it results also completely 
positive, unphysical effects of the kind discussed below would affect this case as well. 
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On the other hand, let us notice that maximal entanglement can be produced without 
any transient. A particularly interesting example is that of two neutral kaons that are 
produced, via the weak interaction, as decay products of a spin-one $ resonance: due to 
angular momentum conservation, the two spin-zero kaons fly apart back to back, in a state 
that resembles that of the singlet for two spin 1/2 particles. [25] 

Using standard techniques, in the next Section we shall derive the Markovian limit 
of the master equation describing two two-level systems in interaction with a stochastic 
environment. After waiting for the correlations in the environment to die out, the resulting 
finite time-evolution F t turns out to be describable in terms of a factorized dynamics: 
= 7t <S>7t, where 7t represent a single open system dynamics, in general non-positive. In 
Section 3 we shall then apply the derived time-evolution F t to a maximally entangled, pure 
initial state, and show that "negative probabilities" may arise even though the dynamics 
7t remains positive on the constituent single-system states. The case of partially entangled 
states and that of mixed entangled states is discussed in Section 4. The concluding Section 
5 contains our final considerations. 



2. MARKOVIAN MASTER EQUATION 

The physical model we shall study is formed by two, non-interacting, two-level systems 
immersed in the same, external environment. The Markovian limit of their subdynamics 
will be derived using the same techniques and approximations widely adopted in analyzing 
single system time-evolutions. [5-9] 

For sake of definiteness, the action of the environment on the two subsystems will be 
assumed to be mediated by a weak time-dependent stochastic field, coupled to their spin- 
like degrees of freedom. [26-28] This choice is of sufficient generality for the considerations 
that follow. Let us point out that this model can describe real physical situations, like the 
ones occurring in interferometric set-ups, involving the propagation of neutrons in random 
magnetic fields [29-36] or photons in random optical media. [5-8, 37, 38] Moreover, it has 
been used to study dissipative effects in correlated neutral mesons under the action of 
weak, stochastic gravitational fields. [39-42] 

Without loss of generality, the total system hamiltonian can be taken to be: [36] 

H = + H<?> + H <" + iff , (2.1) 
B^ = ^*\ H^=±V^(t)*?\ A = 1,2, (2.2) 

2 = 1 

where = &i <g> 1, vf^ = 1 ® Oi, are the two system spin operators, represented by the 
Pauli matrices o % , i = 1, 2, 3, while \^(t) = (v{ A) (t), v£ A) (t),v£ A) (t)), A = 1,2, are the 
stochastic, Gaussian field variables, independently coupled to the spin degrees of freedom 
of the two systems. For simplicity, we assume V( A \t) to have zero mean, (V^ A \t)) = 0, 
and stationary, real, positive-definite covariance matrix [W^ B ^(t)] with entries 

w( AB \t -s) = (v} A) (t)vj B) (s)) = (W^)*(t -s) = W^ A \s - t) . (2.3) 
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For completeness, let us point out that the hamiltonian (2.1) can be equivalently 
interpreted as describing two subsystems in interaction with two independent baths of 
identical physical characteristics; for defineteness, in the following we find more convenient 
to refer to the single bath picture. 

Being coupled to a stochastic field, the complete 4x4 spin density matrix R(t) is 
also stochastic; an effective, "reduced" , spin density matrix p(t) is obtained by averaging 
over the noise: pit) := (R(t)). At the initial time t = we may suppose spin and noise to 
decouple, so that: p = (R(0)) = R(0). 

The dynamical equation for p(t) can be obtained in a standard way from the usual 
Liouville-von Neumann equation for R(t), through the intermediate use of the interaction 
picture. The resulting master equation contains an infinite series of terms. As usually 
done in the case of a single system subdynamics, a simplified, more manageable expression 
for it can be derived by means of physical considerations. [5-9, 26] 

By hypothesis, the action of the external stochastic field on the two subsystems is weak; 
within this "weak coupling limit" assumption, one can then focus on the dominant first 
term in the expansion of the general master equation, neglecting higher order contributions. 
One explicitly finds: 



2 3 



d tP (t) = -i[H^ + ir< 2 \ p(t)\ J2 EC(*)H »H - p{t) \\ ' {2Aa) 

A,B=li,j=l 



3 ,t 



C\f\t) = £ I dsWi AB \s)U kj (s) , (2.46) 
k=i Jo 

where 

/l \ 

Uij(t) — cos u t smu) t (2.5) 




is the orthogonal matrix [Uij(t)] that represents the rotations of the Pauli matrices due to 

the action of the free Hamiltonian: e~ ltH o a^e ltH o } = Y^=i Uij(t) &j A ^- 

Further, by the same physical arguments, the memory effects in (2.4) should not be 
physically relevant: within the abovementioned hypothesis, the use of the Markovian limit 
is therefore justified; in practice, this can be implemented by extending to infinity the 
upper limit of integration in (2.46). More precisely, for situations amenable to a rigorous 
mathematical treatment, one can show that a linear, local in time subdynamics is the 
general result of a limiting procedure in which the coupling constant £ between system and 
external environment, and the ratio r/T between the typical time scale of the system and 
the decay time of the correlations in the environment, become small. [2-4] The quantities £ 
and t/T regulate both the weak coupling limit and the Markovian approximation. 

In order to keep the discussion in the subsequent sections as simple as possible, we 
shall make some further simplifying assumptions on the environment correlations (2.3). 
We first assume that the external stochastic field be oriented along the third direction, 
\( A \t) = (0, 0, V^ A \t)), with exponentially suppressed correlation functions: 

(V 3 {1 \t)V^\s)) = (V 3 {2 \t)vj 2 \s))=g 2 e-^- s \ , (2.6a) 
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(V 3 {1 \t)V^\s)) = fe-^- s \ . (2.66) 

Furthermore, we make the physically sensible hypothesis that the non-diagonal, off-site 
correlations W^ 2 " 1 be subdominant with respect to the diagonal, on-site ones, W^^; in 
practice, this can be achieved by assuming a hierarchy in the strength (f 2 <C g 2 ) and decay 
constants {ji <C v) of the two types of correlations. In this way, the interaction between the 
two subsystems induced by the coupling with the environment becomes negligible. (The 
general case is briefly treated in the Appendix, where more details on the derivation of 
Eq.(2.7) below can also be found.) Then, to lowest order, the finite time-evolution for the 
density matrix, p(0) \— > pit) = r f [p(0)], assumes a factorized form, T t = "ft <S> It, being 
generated by the following Markovian master equation: 

d tP (t)= (L®l + l®L)\p(t)] ; (2.7) 

the linear operator L[-] = L [-] + £i[-], the generator of 7 t , acts on 2 x 2 density matrices 
77, and its explicit form is as follows: 

L o[v] = -i[Ho,rj\ , H = uja 1 , (2.8a) 

L i[v] = «(o"3?7 0-3 -rf)- fl{°~2V 03 + 0-3770-2) , (2.86) 

with 

a = —5 tt , p = —5 tt , u) = — + p . (2.8c) 

In this way, each system evolves independently, with the dynamics generated by^ 

d tV (t) = L[ V (t)] . (2.9) 

This equation is of the Bloch-Redfield type [5-9] and as such it is known not to be positive.* 
As already mentioned in the Introduction, to cure this pathology an ad hoc prescription 
has been proposed, widely adopted in the literature: restrict the possible initial states 7/(0) 
to those for which r)(t) = 7t[//(0)], t > 0, as generated by (2.9), is still a state. As we 
shall see in the next section, this requirement is in general not enough to guarantee the 
consistency of (2.7). 



* It is interesting to notice that essentially the same master equation (2.8), (2.9) is also 
the result of the Markovian approximation of a stochastic dynamical evolution based on 
the " quantum state diffusion" approach. [43, 44] 

* Indeed, for sufficiently small, but positive times, the evolution equation (2.9) will map 

the initial state 77(0) = ( 1 J into a non-positive matrix 77(7;). 
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3. MAXIMALLY ENTANGLED STATES 



As shown in the previous section, the dynamics of two non-interacting systems im- 
mersed in the same bath takes a factorized form, I\ = 7t <S> 7t, at least for times much 
longer than the characteristic correlation times in the environment. However, the initial 
state p(0) of the compound system, on which F t acts, need not be in factorized form: due 
to the short-time interaction with the environment, the subsystems can emerge from the 
transient in an entangled state. 

To avoid inconsistencies with the single-system dynamics 7t, we shall adopt the previ- 
ously mentioned prescription of restricting its action to those states for which positivity is 
guaranteed, for any t > 0. For the compound system under study, this amounts to require 
that the partial traces r/ 1 ) = Tr 2 [/o(0)] and i]^ = Tri[p(0)] over the degrees of freedom of 
the second, respectively the first, subsystem be admissible states for the map 74. 

Let us decompose a generic 2x2 density matrix 77 along the Pauli matrices and the 
identity <jq. 77 = X^=o r i ii(J ^i with rf = 1/2 and rf real; its time-evolution, generated by 
the equation (2.8)-(2.9), is then given by i](t) = 74(77] = Xlu=o 7 7 m (*) c7 M) with components: 



At) 



1 

2 ' 

V 1 (t) = e- 2at V 1 , 



V 2 (t) = e 



-at 



-at 



^ cos 20t - ^ sin 2Q?j rj 2 - sin 2Vtt rf 



^ a 
cos 2m -| — — sin 

—■1 L 



2nt] 7] 3 + i ^- F J- sin 2VLt 7] 2 



n 



(3.1a) 
(3.16) 

(3.1c) 
(3.1d) 



with O = a/w 2 — /3 2 — a 2 /4, and r] 1 , r] 2 , r] 3 the initial density matrix components. 

With the eigenstates r}± = |±)(±| = (cr ± (T\)/2 of the free systems hamiltonian in 
(2.8a), Hq\±) = ±a;|±), one can build the maximally entangled state: 



(3.2) 



Let us assume that after the transient the compound system be in such a state, so that 
the initial density matrix, on which the total Markovian dynamics Y t = ^ t ® It acts, is 
given by: 



P(0) 



V- + V- ® V+ - V+- ® V-+ - V-+ ® V+-) , (3-3) 



with 77 ±T = (0-3 ± ia 2 )/2. The two partial traces = Tr 2 [p(0)] and r/ 2 ) = Tri[p(0)], 
being equal to ctq/2, are left invariant by the dynamics (3.1), and therefore represent 
admissible states for the evolution 7^ 

It is a matter of a simple computation to apply the evolution given in (3.1) to the 
four 2x2 matrices 77+, 77- , 77-1 , 77 |_ and therefore obtain the explicit expression for the 
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evolved 4x4 matrix p(t) = F t [p(0)] = -f t <g> 7t[p(0)]. In the basis for which the Pauli 
matrices assum 
explicitly gets: 



matrices assume the standard form a x = ( ^ 0/ ' 02 = 0/ ,0 " 3= (o 011(1 



Pit) = ~ A 



(A.{t) C(t) C(t) B+(t)\ 

-C(t) A+(t) B_(t) -C(t) 

-C(t) B_(t) A+(t) -C(t) 

\B+(t) C{t) C(t) A_(t)J 



(3.4) 



where 



A±(t) = lie 



-2at 



cos 2Vtt + ^ sin 2Q,t\ + ( - - sin 2 20t 



B±(t) = - e - 4Q *±e- 2at 



cos 2Vtt - ^- sin 20^ + ( - ^-^ ] sin 2 20f 

—■■1 u / \ i L 



C(t) = ie~ 2at sin20t 



cos 2Qt — ^ sin 2Qt 



(3.5o) 
(3.56) 
(3.5c) 



The matrix p{t) in (3.4) should represent the state of the compound system at time 
t, having been originally prepared in the initial entangled state (3.3). The matrix (3.4) 
should then be positive. However, one can easily check that one of its eigenvalues can 
become negative, precisely that corresponding to the eigenvector (1,0, 0,-1). ' Indeed, 
from its expression: 



\(t) = ^ 1 + e~ 4at - 2e~ 2at 



2uj 2 — O 2 
cos 2 2VLt H ^ sin2 2Vtt 



(3.6) 



one checks that A(0) = A(0) = 0, while A(0) = — 8f3 2 , so that X(t) starts assuming negative 
values as soon as t becomes nonzero. 

In order for the map r f = 7 t <S> "ft to produce a physically acceptable dynamics, states 
like (3.3) must therefore be excluded from its domain of definition. Entanglement is crucial 
in revealing this physical inconsistency; indeed, on factorized states r/ 1 ) ®r]^ 2 \ with j/ 1 ), 
1]( 2 ^ admissible starting density matrices for 7 t , the dynamics T t = "ft^lt remains positive. 

The negativity of X(t) will not last for ever: due to the damping factors, the expression 
in (3.6) becomes positive after a certain time, and actually asymptotically tends to 1/4, as 
the remaining three eigenvalues of (3.4). This is a consequence of the dynamics generated 
by (2.7) for which the von Neumann entropy, S[p] = — plnp, always increases (as already 
observed, the map 7*, hence F t , is unital, 7t[co] = cr ); therefore, any initial state p(0) of 



t Using the definitions in (3.5), the four eigenvalues can be explicitly written as: (A+ — 

S_)/4, (A. - B+)/4, {(A+ + A_ + B+ + B_) ± [(A+ -A_+B_- B+) 2 - 16C 2 ] 1/2 }/8. 
The corresponding eigenvectors turn out to be time dependent; their explicit expressions 
are involved and not particularly inspiring. 
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the compound system is asymptotically driven for long times to the maximally disordered 
state p = a <8> cr /4. 

In general, master equations of the form (2.9), (2.86) may involve parameters a and j3, 
not as in (2.8c), but totally independent. In such cases, contrary to (2.8c), a can become 
vanishingly small, without conflicting with the Markovian hypothesis. Consequently, the 
eigenvalue \(t) becomes a periodic function of time and assumes negative values even for 
arbitrary large times. 



4. PARTIALLY AND MIXED ENTANGLED STATES 

As already observed, it is the entanglement of the initial state p(0) of the otherwise 
independent two subsystems that allows revealing the unphysical effect of production of 
"negative probability" by the dynamics IV The magnitude of the phenomenon is directly 
connected to the amount of entanglement that the initial state p(0) contains. 

This can be easily shown by taking the partially entangled state 

\ipo) = cosO |+) <g> |-) - sin0 |-) <g> |+) , (4.1) 

as starting state, instead of the maximally entangled one in (3.2). The evolution in time of 
the corresponding density matrix pg(0) = \ipg)(ipg\ can be easily obtained as before using 
the explicit expressions in (3.1).^ One finds that also in this case the eigenvalue Xe(t) of 
pe(t) = T t [pg(0)} corresponding to the eigenvector (1, 0, 0, —1) can assume negative values; 
in fact, Afl(0) = Xg(0) = 0, while A e (0) = 2(a 2 cos 2 2# - 4(3 2 sin 2 20), which is negative 
provided: tan 2 29 > a 2 /4(3 2 . 

In other terms, once the evolution is given, and therefore the parameters a and (3 
of the corresponding master equation are fixed, the time evolution F t = ^ t <g> ^ t of the com- 
pound system becomes physically inconsistent on initial states that possess a sufficiently 
high degree of entanglement. Therefore, in order for Y t to be an acceptable Markovian 
evolution, one has to further restrict its domain of definition, in order to exclude also those 
partially entangled states. 

The discussion can be extended to entangled mixed initial states, like the Werner 
states: [45] 

1 -p 1 
Pw=PP^ —a o 0a o , ~-<p<l, (4.2) 

where p is again the maximally entangled state in (3.3). Also in this case one can show that 
the eigenvalue \wif) of pw(t) — T t [pw] corresponding to the eigenvector (1, 0, 0, —1) can 
take negative values, provided the parameter p, that measures the degree of entanglement, 
is sufficiently close to one. 

t Notice that as before the partial traces r]^ = Tt 2 [pg(0)] = a /2 + cos26>cri/4 and 

(2) 

i] e = Tri[pe(0)] = o"o/2 — cos26>cri/4 are perfectly admissible states of the dynamics "ft- 
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The discussion becomes particularly transparent when a = 0. In this case, the eigen- 
values of pw(t) become a periodic function of time; for Xw(t) one then explicitly obtains: 



X w (t) = -{l+p 



1-2 cos 2 2Qt + 



2u 2 - n 2 

02 



sin 2 2Qt 



(4.3) 



From this expression, one sees that the minimum value of Xw(t) becomes periodically 
negative provided p > {uj 2 — f3 2 )/(uj 2 + 3/3 2 ). This possibility is excluded in absence of 
noise: (3 = 0; in this case, the evolution 74 (and hence F t = -f t ®7t), being unitary, results 
automatically positive. 



5. DISCUSSION 

The dynamics of a subsystem in weak interaction with an external environment can 
be described in terms of linear maps obeying a Markovian master equation. This result has 
been rigorously proven in very special cases; it is nevertheless believed to hold in general on 
the basis of simple physical considerations: when all correlations in the environment have 
died out, non-linearities and memory effects should disappear from the reduced subsytem 
dynamics. 

As shown in Sect. 2, the same type of arguments allow deriving a Markovian limit for 
the master equation describing the time evolution of two non-interacting subsystems in 
contact with the same reservoir; the corresponding dynamical map F t for the compound 
system turns out to assume a factorized form: F t = ^ft® It- 
Fa. the case of two-level systems, 7$ usually takes a Bloch-Redfield type form, and 
therefore it is not in general positive. To avoid inconsistencies, one usually restricts the 
possible initial states to those for which ^ t remains positive (the so-called "slippage of initial 
conditions" ) . This prescription works also in the case of the evolution F t = ^ t <g> 7 f for two 
identical subsystems, provided the initial state is in separable form: p(0) = J2i Pi Vi ^Vi^ > 
Pi > 0, YliPi = 1) where 77^ and 77P are admissible states for the first and second 
subsystems, respectively. 

On the contrary, as shown in the previous sections, when the initial state p(0) is not 
in factorized form and the degree of entanglement is sufficiently high, the evolved matrix 
pit) = F t [p(0)] fails to be positive for all times. In keeping with the same attitude adopted 
for a single subsystem dynamics 74, to cure this additional inconsistency one can further 
restrict the domain of applicability of 74 <S>7t- However, this is again a temporary solution: 
indeed, the whole discussion needs be repeated when three or more subsystems in contact 
with the same bath are considered: clearly, further restrictions on 7 f need to be imposed. 

These considerations can not be dismissed as being purely academic; on the contrary, 
they seem to have a direct experimental relevance: as mentioned in the Introduction, 
couples of systems in an entangled state are in fact actively studied, and the ongoing 
experiments on correlated neutral kaons constitute a significative example. From this 
perspective, the widely used cure of redefining the initial conditions in case of non-positive 
Markovian dynamics does not appear to be completely satisfactory. 
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In closing, let us mention that in the few cases for which the Markovian limit of the 
subdynamics can be obtained in a rigorous way, the resulting evolution map 7t turns out to 
be not only positive, but also completely positive. [1-4, 11-16] In these cases, the compound 
map r f = 7^ <S> 7t is also completely positive and therefore no inconsistencies can arise, 
even when T t acts on entangled states. [46, 47] 



APPENDIX 



In Sect. 2, we have seen that the master equation that describes the time evolution of 
two, non-interacting subsystems in contact with the same stochastic bath can be written 
in the following closed form: 



d t p(t) = -i\H^+H^ , p(t)} fl E c 



(AB) 



, pit) 



A,B=1 i,j=l 



C. 



(AB) 
ij 



3 „ c 



dsW^ B \s)U kj (s) , 



(Al) 



(A2) 



after the weak coupling limit and the Markovian approximation have been taken into 
account; here, W^ B ^ represent the correlation functions in the environment, while Uij is 
the orthogonal matrix in (2.5) that takes into account the rotation of the Pauli matrices 
generated by the free Hamiltonian. 

From the properties (2.3) of the correlation functions, it follows that the diagonal, 
on-site coefficients C\^ A \ A = 1,2, are real matrices, that can thus be decomposed into 



symmetric, S\f = (C™ + C£ A >)/2, and antisymmetric, = {C^ - C£ A >)/2, 

components. Correspondingly, the second term in (Al), can be decomposed into Hamil- 
tonian and purely dissipative pieces, so that the total master equation can be rewritten 
as: 

d t p(t) = -i[H,p(t)] + L D [p{t)) , (A3) 
where the total Hamiltonian is now given by 



H = + + Hff + H$> , 



H 



(A) 
D 



J2 e vkAfv. 



(A) JA) 
k 



{AA) 



i,j,k=l 



while the dissipative contribution Ld = + + have diagonal and off-diagonal 
pieces: 



ij=l 
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^ = E (c" 2) + C S K 

ij=l 



(A5a) 
(A56) 



Without additional knowledge on the behaviour of the correlation functions W^ B ^ of the 
environment variables, the form of the evolution equation (A3)-(A5) can not be further 
simplified. 

On the other hand, with the assumptions (2.6), 

W^\t -s) = Wif\t -s)= g*e-^-°\ , W& 2 \t - s) = /* e -l*-l , (A6) 

and all remaining entries zero, the form of the coefficients in (A2) can be explicitly com- 
puted, and a more manageable expression for (A3) can be obtained. Indeed, after some 
simple manipulations, one finds: 

d t p{t) = (l <g> 1 + 1 ® L ) [p{t)] + (Li ® 1 + 1 ® Lx) [p(£)] - L 2 [p(£)] , ( A7) 

where Lq ands Li are linear operators acting on 2 x 2 density matrices rj 

L [ V ] = -i [H, r{\ , H=(^ + P)a 1 , (A8a) 
Li[rj\ =a(a 3 r)a 3 - rj) - (3(a 2 r}&3 + ^3 V^z) , (A.8b) 

while L 2 takes the form: 

L>2[p] =7 {°"3 <S> &3 1 p} - cr 3 (g) cr /9 0"o ® cr 3 - cr ® cr 3 /9(J 3 (g) (7 j 

{(T 3 ® cr 2 , p} + {cr 2 <g> 03, p} - cr 3 ® a p cr <8> cr 2 (A8c) 
- cr <g> cr 2 p cr 3 <E> cr - cr 2 <E> <Jq p <J <g) f7 3 - cr <E> cr 3 p cr 2 <g) cr j 



-5 



The four constants a, 7, 5, that measure the relative strength of the various dissipative 
contributions, are determined by the parameters appearing in the correlation functions 
(A6): 

a= 2 ^ 3= g2uo 7= 2pU 6= Pu ° (A 9) 

(12) 

When the strength and decay constants of the off-diagonal correlations W 33 are much 

smaller than the corresponding ones in W 33 A \ f 2 <C g 2 , 1/u <C 1/p, the dissipative 
constants 7 and 5 can be neglected with respect to a and /3, and the resulting evolution 
equation in (A7) reduces to that presented in (2.7) and (2.8). 



11 



REFERENCES 



1. E.B. Davies, Quantum Theory of Open Systems, (Academic Press, New York, 1976) 

2. V. Gorini, A. Frigerio, M. Verri, A. Kossakowski and E.C.G. Surdarshan, Rep. Math. 
Phys. 13 (1978) 149 

3. H. Spohn, Rev. Mod. Phys. 53 (1980) 569 

4. R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Lect. 
Notes Phys. 286, (Springer- Verlag, Berlin, 1987) 

5. W.H. Louisell, Quantum Statistical Properties of Radiation, (Wiley, New York, 1973) 

6. C.W. Gardiner and P. Zoller, Quantum Noise, II ed., (Springer, Berlin, 2000) 

7. M.O. Scully and M.S. Zubairy, Quantum Optics (Cambridge University Press, Cam- 
bridge, 1997) 

8. R.R. Puri, Mathematical Methods of Quantum Optics, (Springer, Berlin, 2001) 

9. H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford 
University Press, Oxford, 2002) 

10. R. Dumcke and H. Spohn, Z. Physik B34 (1979) 419 

11. E.B. Davies, Comm. Math. Phys. 39 (1974) 91 

12. E.B. Davies, Math. Ann. 219 (1976) 147 

13. V. Gorini, A. Kossakowski and E.C.G. Surdarshan, J. Math. Phys. 17 (1976) 821 

14. V. Gorini and A. Kossakowski, J. Math. Phys. 17 (1975) 1298 

15. A. Frigerio and V. Gorini, J. Math. Phys. 17 (1976) 2123 

16. G. Lindblad, Commun. Math. Phys. 48 (1976) 119 

17. A.J. van Wonderen and K. Lendi, J. Stat. Phys. 80 (1995) 273 

18. A. Royer, Phys. Rev. Lett. 77 (1996) 3272 

19. D.A. Lidar, Z. Bihary and K.B. Whaley, Chem Phys. 268 (2001) 35 

20. S. Gnutzmann and F. Haake, Z. Phys. B 10 (1996) 263 

21. A. Suarez, R. Silbey and I. Oppenheim, J. Chem. Phys. 97 (1992) 5101 

22. P. Gaspard and M. Nagaoka, J. Chem. Phys. Ill (1999) 5668 

23. J. Wielkie, J. Chem. Phys. 114 (2001) 7736 

24. L. Diosi, Physica A199 (1993) 517 

25. The Second Daxj)ne Physics Handbook, L. Maiani, G. Pancheri and N. Paver, eds., 
(INFN, Frascati, 1995) 

26. J. Budimir and J.L. Skinner, J. Stat. Phys. 49 (1987) 1029 

27. B.B. Laird, J. Budimir and J.L. Skinner, J. Chem. Phys. 94 (1991) 4391 

28. B.B. Laird and J.L. Skinner, J. Chem. Phys. 94 (1991) 4405 

29. J.L. Staudenmann, S.A. Werner, R. Colella and AW. Overhauser, Phys. Rev. A 21 
(1980) 1419 



12 



30. S.A. Werner and A.G. Klein, Meth. Exp. Phys. A23 (1986) 259 

31. V.F. Sears, Neutron Optics, (Oxford University Press, Oxford, 1989) 

32. Advance in Neutron Optics and Related Research Facilities, M. Utsuro, S. Kawano, 
T. Kawai and A. Kawaguchi, eds., J. Phys. Soc. Jap. 65, Suppl. A, 1996 

33. K.C. Littrell, B.E. Allman and S.A. Werner, Phys. Rev. A 56 (1997) 1767 

34. B.E. Allman, H. Kaiser, S.A. Werner, A.G. Wagh, V.C. Rakhecha and J. Summham- 
mer, Phys. Rev A 56 (1997) 4420 

35. H. Rauch, S. A. Werner, Neutron Interferometry (Oxford University Press, Oxford 
2000) 

36. F. Benatti, R. Floreanini and R. Romano, J. Phys. A 35 (2002) 4955 

37. E. Collett, Polarized Light, (Marcel Dekker, New York, 1993) 

38. C. Brosseau, Fundamentals of Polarized Light, (Wiley, New York, 1998) 

39. F. Benatti and R. Floreanini, Banach Center Publications, 43 (1998) 71 

40. F. Benatti and R. Floreanini, Nucl. Phys. B511 (1998) 550 

41. F. Benatti and R. Floreanini, Ann. of Phys. 273 (1999) 58 

42. F. Benatti, R. Floreanini and R. Romano, Nucl. Phys. B602 (2001) 541 

43. T. Yu, L. Diosi, N. Gisin and W.T. Strunz, Phys. Rev. A 60 (1999) 91 

44. T. Yu, L. Diosi, N. Gisin and W.T. Strunz, Phys. Lett. A265 (2000) 331 

45. R.F. Werner, Phys. Rev. A 40 (1989) 4277 

46. F. Benatti and R. Floreanini, Chaos, Solitons and Fractals 12 (2001) 2631 

47. F. Benatti, R. Floreanini and R. Romano, J. Phys. A 35 (2002) L551 



13 



